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Abstract. We address several problems concerning the geometry of the space of Hermitian 
operators on a finite-dimensional Hilbert space, in particular the geometry of the space of density 
states and canonical group actions on it. For quantum composite systems we discuss and give 
examples of measures of entanglement. 

1. Introduction 

In his book pQ, Dirac uses the description of the interference phenomena, via 
the superposition rule, to justify the requirement of linearity on the carrier space 
of states to deal with quantum evolution. However, the probabilistic interpreta- 
tion of state vectors forces on us the identification of physical (pure) states with 
points of the complex projective space associated with the starting vector space 
of "states". With this identification a (global) linear structure is no more avail- 
able, now interference phenomena will be described with the help of a connection 
(Pancharatnam connection) |2j . The Hermitian structure available on the starting 
Hilbert space of "states" induces a Kahlerian structure on the complex projective 
space. The induced action of the unitary group, projected from the one on the 
Hilbert space, allows for the imbedding of the complex projective space into the 
dual of the Lie algebra of the unitary group itself by means of the momentum map 
associated with the symplectic action of the group. Within this ambient space, by 
means of the available linear structure, it is possible to construct convex combina- 
tions of pure states (rank-one projectors) and build the totality of density states. 
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The space so constructed is not linear and gives rise to interesting geometrical 
structures. To deal with these various non-linearities, recently [3], elaborating on 
previous geometrical approaches to quantum mechanics \4., we have considered 
the differential geometry of density states. This approach seems to be quite ap- 
propriate to deal with composite systems and the set of separable and entangled 
states which do not carry a linear structure. In this note we would like to further 
elaborate on some subtle points which we have encountered in our previous paper 
[3]. To make the paper self-contained we briefly recall the main results from our 
previous treatment. The paper is organized as follows: after introducing notations 
and conventions, in Section 3 we describe the basic geometric structures useful 
in description of the density states, in particular the invariant Kahler structures 
on the orbits of unitary representations as well as an action of the general linear 
group on the dual of the Lie algebra of the unitary group. The description of 
this action in terms of the Kraus operators along with some of their properties 
is further elaborated in Sections 5 and 6, where a general linear group action on 
density states is discovered. In Section 7 we describe the geometry of the set of 
density states as a convex body, in particular we discuss the smoothness of its 
boundary. The rest of the paper is devoted to description of the composite sys- 
tems: in Section 8 we describe the pure and mixed states of such systems in terms 
of the Segre imbedding and give a general prescription for the construction of en- 
tanglement measures. Examples of such constructions are given in Sections 9 and 
10 for bipartite and multipartite systems. 



2. Notations and conventions 

Let 7i be an n-dimensional Hilbert space with the Hermitian product (x, y)^ 
being, by convention, C-linear with respect to y and anti-linear with respect to 
x. The unitary group U(TC) acts on H preserving the Hermitian product and it 
consists of those complex linear operators A £ gliTL) on TL which satisfy AA^ = 
I, where A^ is the Hermitian conjugate of A, i.e., (Ax,y)-H = (x, A^y)^. The 
geometric approach to Quantum Mechanics is based on the realification Hr of 7i 
considered as a Kahler manifold (TCr, J, g, to) with canonical structures: a complex 
structure J, a Riemannian metric g, and a symplectic form uj. The latter come 
from the real and the imaginary parts of the Hermitian product, respectively. 
After the obvious identification of the vectors tangent to Tin with TCn, all these 
structures are constant (do not depend on the actual point of Tin) and read 

J(x) = i-x, g(x,y) +i -u)(x,y) = (x,y) n . 

We have obvious identities 



J 2 = -I, u(x,Jy)=g(x,y), g(Jx, Jy) = g(x,y), u(Jx, Jy) = u(x,y). 
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Fixing an orthonormal basis (e&) of TC allows us to identify the Hermitian 
product (x, y)ji on TC with the canonical Hermitian product on C n of the form 
(a, b)c n = Ylk=i ®k~bk, the group U(TC) of unitary transformations of TC with U(n), 
its Lie algebra u(TC) with u(n), etc. In this picture (ajk)^ = (akj) and (T^T)jk = 
(aj,ak), where = (tjk) £ C n are columns of the matrix T = (tj^). The choice 
of the basis induces (global) coordinates (qk,Pk), k = 1, . . . , n, on Hr by 

(ek,x) n = (qk + i-Pk)(x), 

in which d qk is represented by e k and d Pk by i ■ e&. Hence the complex structure 
reads 

J = ( d P* ® dqk ~ d 1k ® d Pk) , 
k 

the Riemannian tensor 

9 = ^2 ( dqk ® dqk + dpk ® dpk ^ = 2 ( dqk V d # fe + dpfc V dpk ^ 
k k 

and the symplectic form 

u = ^2dq k Adp k , 
k 

where xVy = x(3y + y(&x\s the symmetric, and x /\y = x (& y — y (& x \s the 
wedge product. In complex coordinates Zk = qk + i'Pk one can write the Hermitian 
product as the complex tensor (•, = J^ fc dlk dzk- 

One important convention we want to introduce is that we will identify the 
space of Hermitian operators A = A' with the dual u*(TC) of the (real) Lie algebra 
u(TC), according to the pairing between Hermitian A £ u*(TC) and anti-Hermitian 
T £ u(TC) operators (A,T) = § ■ Tv(AT). The multiplication by i establishes 
further a vector space isomorphism u(TC) 3 T iT £ u*(TC) which identifies the 
adjoint and the coadjoint action of the group U(H), Adc/(T) = UTUl Under this 
isomorphism u*{TL) becomes a Lie algebra with the Lie bracket [-4,5] = \[A, £?]_ 
(where L4,i?]_ = AB — BA is the commutator bracket), equipped additionally 
with the scalar product 

{A, B) u * = ^Tr(AB) (1) 

and an additional algebraic operation, the Jordan product (or bracket) [A, B] + = 
AB + BA. The scalar product is invariant with respect to both: the Lie bracket 
and the Jordan product 



u*(H) — (A [£, 5]) u *Cft)> 
([A i] + , B )u*(H) = (A [€,B] + )u*(H)- 



(2) 
(3) 
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and it identifies once more u*(H) with its dual, u*(7i) 3 A i— ► A = jA S u(H), so 
vectors with covectors. Under this identification the metric Q corresponds to the 
invariant metric 

(A,B) u = ±Tr(AB) (4) 

on u(H) which can be viewed also as a contravariant metric on u*(7i). 

For a (real) smooth function / on Tiji let us denote by gradf and Harrif the 
gradient and the Hamiltonian vector field associated with / and the Riemannian 
and the symplectic tensor, respectively. In other words, <?(•, gradf) = df and 
cj(-, Harrif) = df. We can define also the corresponding Poisson and Riemann- 
Jordan brackets of smooth functions on H by 

{/,/% = Lo(Hamf,Ham f/ ), {f,f'} g = g(grad f , gradf) 

and the total bracket by {f,f'}n = {f>f'}g + Note that any complex 

linear operator A £ gl(T~i) induces a linear vector field A on 7i by A{x) = Ax and 
a quadratic function Ja(x) = ^(x,Ax)n- The function /a is real if and only if A 
is Hermitian, A = A\ It is easy to see the following. 

THEOREM 1. 

(a) For Hermitian A we have 

gradf A = A and Hanif A = i A. 

(b) For all A, B £ gl(7~L) we have {/a, Ib}h = Hab- In particular, 

= Iab+ba, (5) 

{/a,/b}o; = f-i(AB-BA)- ( 6 ) 



The unitary action of U(Tl) on Tt is in particular Hamiltonian and induces a 
momentum map ^ : Hr — ► u*(7i). The fundamental vector field associated with 
jA G uiTL), where A G u*(Ti) is Hermitian, reads i^4, since 

exp ( — A)(x)=iA(x). 
dt\ t=0 i 

The Hamiltonian of the vector field iA is /a, so the momentum map is defined by 

(/j,(x), — A) = f A (x) = -(x, Ax) n . 

But by our convention 

</i(x), X -A) = Z -Tr(»(x)~A) = ^Tv(^x)A), 
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so that Tr(p(x)A) = {x, Ax)^ and finally, in the Dirac notation, 

p{x) =| x)(x\ . (7) 

Hence, the image of this momentum map consists of all non-negative Hermitian 
operators of rank < 1. The space of all non-negatively defined operators, i.e. of 
those p G gl(H) which can be written in the form p = T^T for a certain T G gl(7i), 
we denote by V(H). It is a convex cone in u*(7i). The set of density states V(Tt) 
is distinguished in the cone V{TL) by the equation Tr(p) = 1, so it is a convex body 
in u*{7i). Denote by V k (7i) (resp., T> k {7i)) the set of all non-negative hermitian 
operators (resp., density states) of rank k. In the standard terminology, T> l {TL) 
is the space of pure states, i.e. the set of one-dimensional orthogonal projectors 
| x)(x |, \\x\\ = 1. 

It is known that the set of extreme points of T>(7i) coincides with the set T> l (TL) 
of pure states (see Corollary I2J . Hence every element of T>(7i) is a convex com- 
bination of points from V l (TL). The space V 1 (7i) of pure states can be identified 
with the complex projective space PTL ~ CP™ -1 via the projection 

H\{0}9 I hM^eD 1 (W) 

which identifies the points of the orbits of the C \ {0}-group action by complex 
homoteties. It is well known that P 1 (W) is canonically a Kahler manifold. This 
will be the starting point for the study of geometry of u*(TL) and the set V(TL) of 
all density states. 

3. Geometry of u*(H) 

Recall that u*(TC) is canonically an Euclidean space with the scalar product 
(A,B) U * = iTr(AB). We have also seen that u*(7i) is canonically a Lie and a 
Jordan algebra with the brackets [A, B] = j(AB - BA) and [A, B} + = AB + BA, 
respectively. Note also that, for A being Hermitian, is the pullback = 
p*(A) = A o /i, where A = (A, -) u * = iA G u(7i). The linear functions A generate 
the cotangent bundle T*u*(Tl), so that © and (jHJ) mean that the momentum 
map p relates the contravariant analogs of g and to, respectively, with the linear 
contravariant tensors R and A on u*(TC) corresponding to the Jordan and Lie 
bracket, respectively. The Riemann-Jordan tensor R, defined in the obvious way, 

R(Z)(A,B) = (£, L4, = \Tr(aAB + BA)), (8) 
is symmetric and the tensor 



A(0(A,B) = <£, [A, B]) u * = -Tr(£(AB - BA)), 

1% 



(9) 
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is the canonical Kostant-Kirillov-Souriau Poisson tensor on u*(7i). They form 
together the complex tensor 

(R + i- A)(0(A B) = 2(C, AB) U , = Ti^AB) (10) 

and the momentum map relates this tensor with the Hermitian product. 

On u*{Ti) consider the (generalized) distributions D\ and Dr induced by the 
tensor fields A and R, respectively. To be more precise, Denote by J and TZ 
the (1, l)-tensors on u*(H), viewed as a vector bundle morphism induced by the 
contravariant tensors A and R, respectively, J, R : Tu*(H) — ► Tu*(H), where 

Jj{A) = [A,S]=A € (A), 
TZ^A) = [A,£] + = Rf(A), 

for A £ u*{H) ~ Tgu*(H). The image of J is D A and the image of TZ is D R . It is 
easy to see that the tensors J and TZ commute and 

J^o7Z^A)=7Z^oJ^A) = [A,f]. (11) 

We will consider also the generalized distributions D\ = Dr + D\ and Dq = 
D R f]D A . 

1. Distribution Da- 

It is well known that the distribution D\ can be integrated to a generalized foliation 
T\ whose leaves are orbits of the canonical U (7Y)-action on u*(H) given by U(H) x 
u*(Ti) B (U,£) ' — ^ U^W G u*(H). They are represented by the spectrum of the 
operator, i.e. p and p' belong to the same [/(W)-orbit if and only if they have 
the same set of eigenvalues (counted with multiplicities). Moreover, the orbits 
are symplectic leaves of the Poisson structure A with the corresponding U(H)- 
invariant symplectic form 7]°. These symplectic structures can be extended to 
canonical Kahler structures as shows the following. 

THEOREM 2. 

(a) J"^ is a self-adjoint with respect to (•, ■)„* and non-positively defined operator 
on u*(Ti) with the kernel Da(£,) ± - 

(b) The (1, 1) -tensor J on u*(7i) defined by 

j ifAeD^H) 1 - 
JdA) = \J^ ( A ) V A e (12) 

satisfies = — J and induces an U{TL) -invariant complex structure J on 
every U(H)- orbit O. 
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(c) For every U(Tt)-orbit O the tensor 

1 f(A,B)= v f(A,J^(B)) (13) 

is an U '(H) -invariant Riemannian metric on O and 

1 °(Jt(A),B) = r l f(A,B), A,BeD A (0. (14) 

In particular, (O, J, rf , 7 ) is a homogeneous Kdhler manifold. Moreover, 
if£& u*(H) is a projector, £ 2 = £, and (GO, then = and "f°(A, B) = 
(A,B) U *. 

Remark. The tensor J is canonically and globally denned. It is however not 
smooth as a tensor field on u*(H). It is smooth on the open-dense subset of 
regular elements and, of course, on every U (H)-orbit separately. 

2. Distribution Dr. 

We have some similar results for the tensor 1Z which however are not completely 
analogous, since the distribution Dr is not integrable. 

THEOREM 3. 

(a) 7£| is a self-adjoint with respect to (-, -) u * and non-negatively defined operator 

on u*(H) with the kernel Da(C)" L - 

(b) The (1, 1) -tensor TZ on u*(H) defined by 

v(A\-i^ ifAeDnti) 1 - 

'^ {A) \n^o\(n^ DR ^\-\A) ifAeD R (o 

satisfies TZ 3 = TZ. 



3. Distribution Dq. 

The distribution Dq = D\ f) Dr can be described also as the image of o 7Z% = 
1Zz o In other words, D (0 = {[A,£ 2 ] : A G u*(H)}. 

THEOREM 4. The distribution Do is integrable and the corresponding foliation 
is U(H) -invariant, J -invariant and TZ-invariant, so that J and 7Z induce on 
leaves of T§ a complex and a product structure, respectively. The leaves of the 
foliation T§ are also canonically symplectic manifolds with symplectic structures 
being restrictions of symplectic structures on the leaves of Fa, so the leaves of T§ 
are Kdhler submanifolds of the U(H)-orbits in u*(H). 
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Note however, that Dq coincides with D\ on V(Tt), so that on density states the 
leaves of J-q are just the orbits of the unitary group action. 

4- Distribution D\. 

The distribution D\ = D\ + Dr is the largest one carrying the most qualitative 
information. It turns out to be related to a GL(W)-action as shows the following. 

THEOREM 5. 

(a) The generalized distributions D\ is involutive and can be integrated to gener- 
alized foliations J-\ whose leaves are the orbits of the GL(7i)-action GL(7i) x 
u*(H) 3 (T,f) i-> T£Tt G u*(H). 

(b) The Hermitian operators p and p 1 belong to the same GL(TC)-orbit if and 
only if they have the same number k + of positive and the same number k- 
of negative eigenvalues (counted with multiplicities). Such an orbit, denoted 
by u k+ k (Tt), is of (real) dimension 2nk — k 2 , where k = fc+ + k^, and its 
tangent spaces are described by the formula 

B G ?K +)fc _ Vx,y G KeriO [(Bx,y) n = 0]. (16) 

(c) Any GL(TC)-orbit intersecting V(7i) lies entirely in V(7i), so that V(7i.) is 
stratified by the GL(7i)-orbits. The GL(TC)- orbits in ViTL) are determined 
by the rank of an operator, i.e. they are exactly V k {TL), k = 0, 1, . . . ,n. 

The proofs of all results in this section can be found in 



4. Explicit coordinates on GL(H)-orhits 

A choice of an orthonormal basis in H gives an identification u*(7i) ~ u*(n). 
Let J = {ji, . . . , jk} C {1, . . . , n}. Denote u*j(n) - the open subset in the set u* k (n) 
of rank- A; Hermitian matrices consisting of matrices (a^) £ u k (n) for which the 
submatrix (a rs ) rjSe j is invertible with the inverse (a rs ) rjSe j. 

LEMMA 6. Let £ £ u*j(n). Then the matrix £ is uniquely determined by its rows 
indexed by J, according to the formula 

E rs 

r,s£j 

THEOREM 7. The maps 

$J : u*j{n) -> R k X C^nk~k^-k)l2 _ K 2nk~k 2 ^ 
^^(( a ii)^=l) = (( a n)jgJ; (ars)r<s,sej) 



form a coordinate system in the manifold u k (n). 
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Example. For n = 3 and J = {1,2} 

d\ a\ — ia 2 b\ — ib 2 
a\+ia 2 d 2 c\ — ic 2 ) = (d\, d 2 , a\, a 2 ) G R 4 
b\ + ib 2 ci + ic 2 d 3 

According to the above lemma, the coordinates (c?i, d 2 , a\, a 2 ) are sufficient to 
determine the Hermitian matrix of rank 2 whose principal minor is non- vanishing. 



5. Kraus operators 

The above defined GL("H)-action on u*(TC) is an action by means of particular 
Kraus operators which, by definition, are operations Ka of the form p ^ ^Li AipAj 
for certain A = (Aj) 6 gl(TC) m . It is well known that, considered as operations on 
gl(7i), they are exactly the linear operations preserving the Hermicity and exactly 
the linear operations which are completely positive (cf. [5]). Note that we do not 
assume any normalization condition for Ka, like ^j^l^i = I- It * s eas Y to see 
that the composition of Kraus operators is again a Kraus operator, so that we 
can speak of the semigroup of Kraus operators. Indeed, Ka ° Kb = Ka b, where 
K A . B (p)=Z i , j (A i B j )p(A i B j y. 

Let us note that the space gl(TC) of all complex linear operators on 7i is canon- 
ically a Hilbert space with the Hermitian product (A,B) g i = tjTt(A^B). Using 
the Jamiolkowski isomorphism one can identify the Kraus operator Ka with a 
Hermitian operator on the Hilbert space gl(TC), defined as Pa = J2%PAi, where 
PAi = |Aj) (Aj | . Using the spectral decomposition one can easily see that any Kraus 
operator Ka can be written in a canonical form Ka{p) = Kc{p) = Y^k^pC ■> 
where the operators Ck £ gl(J~L) are pairwise orthogonal, (Cfc, Ck>) g i = for k ^ k', 
and that in this case Aj = aikCk with J2 i aikcqjj = 5^ . 

It is interesting that the operators of the GL("H)-action form exactly the largest 
subgroup in the semigroup of Kraus operators. Since in the literature we could 
find the analogous fact only for the unitary group and the semigroup of normalized 
Kraus operators, we will give a short proof of this general result. 

THEOREM 8. If a Kraus operator Ka is invertible inside Kraus operators, then 
K A {p) = ApA^ for certain A £ GL(H). 

Proof. Assume K^ 1 = K&, so that Yli j(AiBj)p(AiBj)^ = p. According to the 
last observation, AiBj = ctijl for certain aij G C with Ylij l a «j| 2 = 1- There is 
a *oio 0, so Aj and Bj are invertible. Moreover Aj is proportional to -B" 1 , 
namely Aj = a.i j Bj\ Since we can clearly assume that all Aj and all Bj are 
non-zero, we get that all Aj and all Bj are invertible. Indeed, Aj is not invertible 
implies that AiBj = ctij I is not invertible thus zero, so Aj = 0, because Bj is 
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invertible, and similarly for Bj. Moreover, every A{ is proportional to -B" 1 , whence 
to A io , Ai = ji A io , 7i / 0. We get therefore 

k a ( p ) = J2 A iP A l = (ViK)p(ViK)^ 

i 

where 7 = ]T\ |7j| 2 . I 



6. Kraus and GL(H)-action on density states 

We know already that V k (H) are GL(W)-orbits 

in u*(Tt), so smooth submani- 
folds. We used this fact in [3| to show that their intersections with the hyperplane 
Tr(p) = 1, i.e. D k {TL) are smooth submanifolds too. Here we want to stress that, 
m fact, D k (H) can be regarded again as GL(J-L)-oih\ts with respect to an action 
of GLiTi) on T>(TC). Of course, we cannot apply directly the action on u*(7i), 
as T>(Ti) is not an invariant set under this action. We can, however, modify the 
GL(W)-action (in fact, even the Kraus action) on ViTL) in such a way that we get 
an action on density states whose orbits are D k (H). 

For, suppose K A (p) = Y^i A iP A \ is a non- degenerate Kraus operation, i.e. 
Ei A l A i e GL(H). Now, let us define an operation K A on T>(TC) by 

km- KaM 



1r(K A (p)Y 



The definition makes sense, since Tr(K A (p)) > for any density state p. Indeed, 
since 

Tt(K a ( P )) = Tr(J2 ApA\) = Tr((£ A l A i)p) 

i i 

and T = £\ A \ A i is invertible Hermitian and non-negative, so strictly positive, 
and since p in Hermitian non-negative, Tr(Tp) < only if p = 0. 

Now, it is a fundamental observation that we get in this way really an action 
of the semigroup of Kraus operations, i.e. that K A o = K AB . But it is 
straightforward, since 

K B ( P ) 

Tr(K A («^)) Tr(K A oK B (p)) 



7> % ( , A ^b(o)) K A oK B (p) ~ 

Ka ° K-b(p) = x \ = — — = Aa.b(p). 



Tv(K B (p)) 

Note that, though this action is not affine, convex sets are mapped into convex 
sets, as 

K A (\p + (1 - A)p') = \K A (p) + (1 - X)K A (p'), (17) 
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where _ 

T_ _ ATr(if A (p)) 

XTv(K A (p)) + (l-X)Tr(K A (//))' 

Of course, we can reduce this action of the semigroup of non-degenerate Kraus 
operations to its largest subgroup, i.e. to GL(7i), obtaining the action 

GL{H) x V{H) 3 (A,p) ^ A{p) = ^^ f) E V(H). (18) 

This action preserves the rank and one can easily derive from Theorem [T] the 
following. 

THEOREM 9. The decomposition of the convex body of density states T>(TC) into 
orbits of the GL(TL)- action is exactly the stratification V(7i) = \J^ =1 T> (TL) 
into states of a given rank. 

7. The geometry of density states 

The boundary dD(TC) of the convex body of density states consists of the 
states of rank< n, &D(H) = \Jl~l{V k (H). Each stratum V k (H) is a smooth 
submanifold in u*(7i). However, the boundary dT>(TC) is not smooth (except for 
the case n = 2), since its maximal stratum V n ~ 1 (H.) is sewed up along (Jfe=i T> k {7i) 
with edges there, as shows the following (|5J, Theorem 2). 

THEOREM 10. Every smooth curve 7 : R — > u*(TC) through the convex body of 
density states is at every point tangent to the stratum to which it actually belongs, 
i.e. 7 (t) E V k (H) implies T 7 (t) E T l{t) V k {H). 

The above theorem means that inside T>(TC) we cannot smoothly cross the stratum 
T> k (7i) of the boundary transversally to it, like living on a cube we cannot smoothly 
cross an edge of the cube tranversally to it. 

Corollary 1 The boundary of the convex body T>(7i) of density states is a smooth 
submanifold of u*(7i) if and only if dim{7i) < 2. 

Remark. It is well known that for n = 2 the convex set of density states is the 
three-dimensional ball and its boundary - the two-dimensional sphere (so called 
Bloch sphere), so it is a smooth manifold. 

The next problem concerning the geometry of density states we will consider is the 
question of the faces of V{TL), i.e. the intersections of T>(H) with supporting affine 
hyperplanes. In other words, a non-empty closed convex subset Kq of a closed 
convex set K is called a face (or extremal subset) of K if any closed segment in 
K with an interior point in Kq lies entirely in Kq; a point x is called an extremal 
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point of K if the set {x} is a face of K. For p E D(H) consider the decomposition 



W = ?4©Wg = Im(p)0Ker(p), x = < + x p , (19) 
into the kernel and the image of p. 

THEOREM 11. The face ofV(H) through p E V k (H) consists of operators A E 
T>(H) which, according to the decomposition \ltJ\) . have the form A{x p + + Xq) = 
A p + {x p + ) for certain A p + E V(Im(p)), so it is affinely equivalent to the convex body 
of density states in dimension k. 

Corollary 2 Extremal points ofT>(H) are exactly pure states. 

Corollary 3 All non-trivial faces ofT>(H) of maximal dimension, i.e. faces through 

T> n ~ l (H), are tangent to the sphere S(I/n;r), centered at I/n with the radius 

r = , 1 , at points which are collinear with the center and one of the pure 
^/n(n-l) 

states. 

8. Composite systems and separability 

Suppose now that our Hilbert space has a fixed decomposition into the tensor 
product of two Hilbert spaces H = TL 1 ®Ti 2 (of dimensions n\ and ri2, respectively). 
This additional input is crucial in studying composite quantum systems and it has 
a great impact on the geometrical structures we have considered. The rest of this 
paper will be devoted to related problems. 

Observe first that the tensor product map 

(g) -.n 1 xh 2 -^n = n 1 ®n 2 (20) 

associates the product of rays with a ray, so it induces a canonical imbedding on 
the level of complex projective spaces 

Seg : PH 1 x PH 2 -»■ PH = P{H l ® H 2 ), (21) 
{\x l ){x\\x 2 ){x 2 \) ^ \x l ®x 2 ){x l ®x 2 \ . (22) 

This imbedding of product of complex projective spaces into the projective space of 
the tensor product is called in the literature the Segre imbedding [5] . The elements 
of the image Seg(PW 1 x PH 2 ) in PH = PiH 1 (8) H 2 ) are called separable pure 
states (with respect to the decomposition H = H 1 ® H 2 ). 

The Segre imbedding is related to the (external) tensor product of the basic 
representations of the unitary groups UiTL 1 ) and U(H 2 ), i.e. with the representa- 
tion of the direct product group in H = H 1 (8> H 2 , 

UiH 1 ) x U(H 2 ) 3 (p\p 2 ) i ^ p 1 p 2 E U(H) = U(H l H 2 ), 
(p 1 <g> ^(x 1 ® x 2 ) = p 1 {x l ) ® p 2 (x 2 ). 
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Note that p 1 (g) p 2 is unitary, since the Hermitian product in TL is related to the 
Hermitian products in TL 1 and TL 2 by 

{x 1 ®x 2 ,y 1 ®y 2 ) n = (x\y l ) H i ■ {x 2 ,y 2 ) H 2. (23) 

The above group imbedding gives rise to the corresponding imbedding of Lie alge- 
bras or, by our identification, of their duals, which, with some abuse of notation, 
we will denote by 

Segru*^ 1 ) x u*(H 2 ) ^ u*(H), ~ £l ® &■ (24) 

The original Segre imbedding is just the latter map reduced to pure states. In 
fact, a more general result holds true. 

Proposition 1 The imbedding maps V k (7i 1 ) x V l (H 2 ) into V kl (H). 

Let us denote the image Seg(2? 1 (Ti 1 ) X P 1 (W 2 )) - the set of separable pure states 
- by S 1 (Tl), and its convex hull conv (5 1 (W)) - the set of all separable states in 
u*(H) - by S(H). The states from 

8(H) =V{H)\S{H), 

i.e. those which are not separable, are called entangled states. It is well known (see 
e.g. that S 1 (H) is exactly the set of extremal points of S(Tt). 

Of course, by means of the tensor product of representations the group product 
GL(H 1 ) x GL(H 2 ) is canonically embedded in GL(H) like in the case of the unitary 
groups. The canonical actions l|18jl on P(W 1 ) and T>(TL 2 ) give rise to the action to 
the corresponding action of GL(W 1 ) x GL{H 2 ) on V(H l ) x V(H 2 ): 

(A 1 ,A 2 )( Pl ,p 2 ) = (A 1 ( Pl ),Mp 2 )). 
On the other hand, G^H 1 ) x GL(H 2 ) as being embedded in GL(H) acts on V(H). 

THEOREM 12. The aforementioned actions ofGL(H 1 ) x GL(H 2 ) are equivariant 
with respect to the Segre map 

Seg : ViH 1 ) x V{H 2 ) -> V(H). 

Moreover, the set S 1 (Tl) of pure separable states and the set S(7i) of all separable 
states are invariant with respect to the canonical GLiTL 1 ) x GLiTL 2 )- action on 
V{H): 

. ... (A-i (» A->)o(Ai (X) Ai)^ 
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Proof. We get easily 



(A 1 ,A 2 )( Pl p 2 ) 



(A lPl A\) ® (A 2 p 2 Al) 



(A lPl A\) <8> (A 2 p 2 Al) 



TV {[A lPl A\) ® (4 2 p 2 4)) 
3i(pi)(8)^(p 2 ) 



TV(^ 1/0l 4)TV(^ 2/02 4) 



that proves equivariance and the invariance of S 1 (TL). The invariance of S(TL) is 
not automatic, like in the case of the group UiTL 1 ) x U(TC 2 ), since the GL(TC 1 ) x 
GL(7^ 2 )-action is not affine. On the other hand, (j!7j) implies that the convex hull 
thus separability is respected: 



{A 1 ,A 2 ){\pi® P2 + X'p'i® p' 2 ) = \{A 1 ,A 2 ){pi ® p 2 ) + X'(A 1 ,A 2 )(p 1 (g> /9 2 ) 

= Ali(pi) ®i2(p 2 ) + A'2i (pi) (// 2 ), 



where A' = 1 — A and A' = 1 — A. I 

It is a common opinion that the U(H 1 ) x [/(:ft 2 )-action, as preserving separability, 
is crucial for understanding the entanglement. We see, however, that the GLiTL 1 ) x 
GL(W 2 )-action preserves the separability as well and, since the orbits are larger, 
it carries more qualitative information, thus information which is easier to handle. 
As an example, consider the corresponding U{TL l ) x U{TL 2 ) and GL(TL 1 ) x GL(TL 2 ) 
orbits inside pure states. 

Proposition 2 With respect to the Schmidt decomposition 



of the unit vector |W) G TL representing a pure state, the U(Ti ) x UiTL 2 ) -orbits 
in V 1 (TC) are distinguished by the sequence Ai > . . . A m > ; while the GL{T0-) x 
GLiTL 2 ) -orbits are distinguished only by the Schmidt number m. 

Proof. - It is clear that the form of the Schmidt decomposition is preserved by the 
U(TL 1 ) x U(TC 2 )- action. On the other hand, if we have two such decompositions 
Y^k=i ^k\f\) ® \ fl) and XX=i ^k\vl) ® I 7 ?!)' t nen ! srnce fl are pairwise orthonor- 
mal, since <p\ are pairwise orthonormal, etc., there are U l G U(TL l ), i=l,2, such 
that U x (ip\) = r}\ and U 2 {<p 2 k ) = n\, k = 1, . . . , m, so 



m 




(25) 



fc=i 



m m 



(U 1 U 2 )(Y^ hU) ® Wl)) = £ A fc |^) ® Iry 2 ). 



fc=l k=l 



A similar reasoning gives the description of GLiTL 1 ) x GL(W 2 )-orbits, if we only 
take into account that the exact values of the coefficients A& are irrelevant for this 
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action, as the group does not respect the length of the vector. I 

The entangled states play an important role in quantum computing and one of 
main problems is to decide effectively whether a given composite state is entangled 
or not. An abstract measurement of entanglement can be based on the following 
observation [3] (see also Ref. [7j). 

Let E be the set of all extreme points of a compact convex set K in a finite- 
dimensional real vector space V and let Eq be a compact subset of E with the 
convex hull Kq = conv(i^o) C K. For every non-negative function / : E — » R+ 
define its extension fx ■ K — > R + by 

f K (x)= inf (26) 

where the infimum is taken with respect to all expressions of x in the form of 
convex combinations of points from E. Recall that that, according to Krein- 
Milman theorem, K is the convex hull of its extreme points. 

THEOREM 13. For every non-negative continuous function f : E — > R + which 
vanishes exactly on Eq the function fx is convex on K and vanishes exactly on 
K 

Corollary 4 Let F : T> 1 (TL 1 (3TL 2 ) — ► R + be a continuous function which vanishes 
exactly on S l (H l &> H 2 ). Then 

V = *2>(«W) = ViH 1 n 2 ) -» R+ 

is a measure of entanglement, i.e. fi is convex and /i(x) = 44> x £ S(7i). More- 
over, if f is taken U(Ti 1 ) x U (Tt 2 ) -invariant (resp. GLiTi 1 ) x GL{TL 2 ) -invariant), 
then n is [/(H 1 ) x Li \H 2 ) -invariant (resp. GL(H 1 ) x GL(H 2 )) -invariant. 

9. Bipartite entanglement 

One of measures constructed according to the above described prescription is 
the concurrence introduced originally as an auxiliary quantity, used to calculate 
so called entanglement of formation of 2 x 2 systems |HJ. For pure states (|23)l it is 
defined as 

c (*) := y/l-Tr^ = ^||*|| 2 - ^(^(1^1) • Tr 2 (|*><*|)), (27) 

where TVj, i = 1,2, denotes tracing over the i-th. subsystem, and p\ := Tr2 1 1 3 / )( lL I f | . 
It is clear that, indeed, it vanishes for separable states (for which r Tr 1 p 2 = 1). For 
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further generalizations (see Section 10) it is convenient to rewrite l|27[) in slightly 
different form. To this end let us define: 

A:H®H ^H®H, A = 4Pl 1} ® P {2] , (28) 

(i) 

where P_ is the orthogonal projection on the antisymmetric part 7i l A 7i l of the 
tensor product TC l <8> ?T, and we identify in an obvious manner Tt 1 (g) TL 2 <g> Ti. 1 (g) TL 2 
and TL 1 <g> TL 1 <8> W 2 <g> W 2 . It is now a matter of a straightforward calculation that 
c(^) can be expressed as: 



c (¥) = yJ((V\ <g> (tf|)A(|tf> <g> (29) 
Extension of the concurrence l|27jl to mixed states is defined via (|26|). i.e. 

c(p)= inf r Vt,c(4 (30) 

Calculation of c(p) for an arbitrary mixed state p requires a high dimensional 
optimization procedure, it is however possible to derive a lower bound for it which 
in general suffices to discriminate between a separable and an entangled state. 
Our bound is given by a purely algebraic expression easily evaluated for arbitrary 
states and can be, if needed, tightened numerically by optimizing over a parameter 
space of much lower dimensionality than it is demanded by the original definition 
(|30|) [H1E3]- To this end we first replace the by the subnormalized states 
\ipi) = y/~ti\*&i) in Eq. (|30|) . Given a valid decomposition p = ^ \4>i)(cf)i\ into 
M subnormalized states {\(f>i),i = 1,...,M}, any other suitable set {\ipi),i = 
1, . . . , N} such that 

P = X>iXV<i|, (31) 

is obtained by: 

M N 
3=1 i=l 

where both N and M are not smaller than the rank r of p. It can be shown ^2] 
that for the purposes of the present considerations, it is enough to take N < n\n 2 - 
We can now choose as the starting point eg. the decomposition {\4>i), % = 1, . . . , r} 
of p defined in terms of its (subnormalized) eigenvectors 

r 

P = ^2\4>i)(4>i\, p\<t>i) = <H\<t>i), (33) 
i=i 
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where a,i,i = l,...,r, are non-vanishing eigenvalues of p. Now the concurrence 
can be rewritten as: 



c( P ) 



i^^/[(V0V)A(V^V^]: 



where 



Aim 



((&|®(&»l)^(l0i>®l&», 



and the infimum is now taken on over matrices V fulfilling (§32 
simplifies further if expressed in terms of eigenvectors of A, i.e. A 
a = 1, . . . , m = n\(ni — 1)^2(^2 — Namely 



c( P )=mfJ2JE\l VTaV ^i 



(34) 
(35) 

The expression 

Z;ct=l IXa)(Xa| > 

(36) 



where T" k 



(XaKl&j) ® \<f>k})- Obviously any given decomposition (|31l) provides a 



straightforward upper bound of the concurrence, c(p) < Y2i y l^i? (VOI • From 
the point of view of distinguishing separable and entangled states it is much more 
interesting to find a lower bound for c in an an effective way. To this end let us 
write, using Cauchy-Schwarz inequality, 



JE I ^ TayT ^ I JE hi 2 ^ E I [v^y^u I > I E 

y a y a a a 

for arbitrary z a , a = 1, . . . , m. We obtain thus: 



c(p) > inf 
v 



N 

E|[ 



V 



(37) 



(38) 



for arbitrary z a such that ^ Q |z a | 2 = 1. The infimum over V can be effectively 
performed and is given by max {Ai — ^>i ^ 0}> where Xj are the singular values 
of T = Yl a z aT a , i.e. the square roots of the eigenvalues of the positive hermitian 
matrix TT^ in the decreasing order (T^j. The obtained bound still depends on 
the choice of the z a , what allows to tighten the estimate. Thus, one is left with 
an optimization problem on an 2m-dimensional sphere. Note that the constraint 
I -2a 1 2 = 1 is by far simpler to implement than V'V (cf. I32j) . Moreover, the 
dimension m of optimization space is significantly reduced as compared to the 
dimension n\n^ of the original optimization problem defined by Eq. I|32j) . Let us, 
however, point out that any choice of z a gives some lower bound and taking eg. 
all but one z a equal to zero, we can dispose of the optimization entirely if we are 
only interested whether c is positive, which is enough to establish nonseper ability 
of the state in question |Hj . 
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10. Multipartite entanglement 

Separability of multipartite systems, where the Hilbert space of the whole 
system Ti has a fixed decomposition into the tensor product of Hilbert spaces 
7i = Ti 1 (g> . . . , 7i K of subsystems of dimensions m, . . . , Uk is defined by a straight- 
forward extension of the two-component case (cf. Section 8), i.e. via a canonical 
imbedding of the product of projective spaces PTi 1 x . . . x P7i K into the projectivi- 
sation of the tensor product P{TL l (8> • • • TL K ) and the corresponding imbedding 
on the level of Lie algebras and their duals. The pure separable states are thus 
identified with the image under this imbedding of V 1 (TL 1 ) x . . . x V 1 {H K ) and its 
convex hull with the set of all separable states. 

In order to investigate the separability of multipartite states we proposed the 
following generalization of the concurrence considered in the preceding section |14j . 
Let us, in an analogy with ()28j) define 

K 

A {s .y.H®H^n®U, A {s . } = 2 K (g)pM, (39) 

i=i 

where Sj = ± and P^ (respectively p£ '"') are orthogonal projections on the anti- 
symmetric (resp. symmetric) subspace W AW (resp. W J VW J ), define concurrence 
for pure states as 

c {sj} (*) = J((*\®(*\)A {8j} (\*)®\*)), (40) 
and its extension to mixed states by 

C ^= P ^^ UC{ ^ m - (41) 
Closer examination of the action of A{ s j. reveals that if an odd number of pro- 
jectors on anti-symmetric subspaces appears in its definition, the corresponding 
c {sj} 0^0 vanishes identically. Moreover, if S{ = + for all i, i.e. when only pro- 
jections on symmetric subspaces are involved, A^ Sj y is not helpful in detecting 
entanglement [141 115j . 

The techniques which were devised to ease the task of estimating the con- 
currence for arbitrary states in the bipartite case in the previous section, can be 
generalized in a straightforward manner, because the algebraic structure of the 
above ET-partite concurrences is strictly identical to the bipartite definition ([30(1 . 
Thus one can invoke the Cauchy-Schwarz and the triangle inequality and bound 
the concurrence of an arbitrary mixed state from below by 

i V a i 

= maxjAi-^Aj.o}, (42) 
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where T = Y. a z aT a , Tf k = (x a \(\<f>j) ® |0fc}), A {s . } = ^ Q=1 \xa)(Xa\, \<j>k) ar e 
the subnormalized eigenvectors of p, V defines the transition between different 
decompositions of p as in (|32|) . and Xj are eigenvalues of T^T in decreasing order. 
For the proof of the last equality in (|42j) see or |1U| . As in the bipartite case, 
the inequality (|42|) holds for an arbitrary set of complex numbers z a , such that 
S« \ z a\ 2 = 1) what allows for further optimization. 

The above does not only apply to the discrete set of concurrences discussed so 
far, but also to the following continuous interpolation between them: Instead of a 
single direct product of projectors onto symmetric and antisymmetric subspaces, 
one may equally well consider convex combinations thereof, 

A = 2 K £ PS1 _ SK PW®...®PW (43) 

Sl,...,S;v 

where Sj G {+,—}, Psi,...,s x > and the summation is restricted to contributions 
with an even, non-zero number of projectors onto anti-symmetric subspaces. The 
corresponding pure-state concurrence ()40|) can be written in terms of the partial 
traces: 

c(*) = / «sTr(Tr A (|*)<*|)2), (44) 

V se2{ 1 >-- K } 

where 2^- 1 '—' K ' denotes the set of all subsets of {1, ... ,K}, and 

as = P«i.-.«jc II Si ' ( 45 ) 

S!,...,s K igS 1 

Various choices of the coefficients p Slv .. iSx allow to distinguish different categories 
of multipartite entanglement. As an illustration, let us focus on some exemplary 
tri- and four-partite concurrences. 

The so called biseparable pure states (i.e. states taking the form of tensor 
product of a state of one subsystem with a, possibly entangled, state of the other 
two subsystems) in the tri-partite case are easily detected with A = P^"' tg) P^ (g) 
pi 3) , A = <g> pj 2) <g) p| 3) , and A = P^ <g> p| 2) <g> pj 3) . Whereas corresponding 
concurrences vanish identically for bi-separable states like = <8> IC3) f° r 
the first and second choice of A, the last one which reduces to the bi-partite 
concurrence of | dis- 
similarly, different kinds of separability are also captured in larger systems. 
For example, concurrences defined with the help of A = 4Pj^ <g> pjj (g) p] 3 -* (g) Psf , 
Si = Sj = +, and = — for i 7^ k 7^ j, determine with respect to which bipartite 
partition a mixed 4-particle state is separable, and quantify the value of bi-, resp. 
tri-partite concurrences of the entangled part |14j . 
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